Theory of electromechanical coupling in dynamical graphene 
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We study the coupling between mechanical motion and Dirac electrons in a dynamical sheet of graphene. We 
show that this coupling can be understood in terms of an effective gauge field acting on the electrons, which has 
two contributions: quasistatic and purely dynamic of the Berry-phase origin. As is well known, the static gauge 
potential is odd in the K and K' valley index, while we find the dynamic coupling to be even. In particular, the 
mechanical fluctuations can thus mediate an indirect coupling between charge and valley degrees of freedom. 
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Graphene has become one of the most remarkable existing 
material: While the underlying theory is nonrelativistic, its 
electronic quasiparticles appear to exhibit relativistic features 
as a consequence of the low-energy Dirac spectrum lUmi. 
This makes the emergent Dirac fermions Lorentz rather than 
Galilean invariant f6 |. However, this is extremely fragile, as 
any coupling beyond the eflFective model breaks this Lorentz 
symmetry, as does, for example, Coulombic interaction be- 
tween electrons or scattering of Dirac electrons by localized 
impurities. Indeed, a more complete theory must, to a good 
approximation, be Galilean invariant, which is made explicit 
if we consider not only the electrons in the frame of reference 
of a perfect static lattice but the entire graphene membrane as 
a combined dynamical system. 

One prominent feature of the dynamical graphene is lat- 
tice dragging [71, which could be understood in terms of the 
Berry phases that electrons accumulate when moving through 
a time-dependent lattice. The electron dragging due to lattice 
translations depends on the difference between the bare elec- 
tron mass m and the effective electron mass in the crystal m* 
Ul. The larger the difference the stronger the dragging, such 
that these effects should ceitainly be relevant in graphene, 
where the effective electron mass vanishes flSJ. More gener- 
ally, such BeiTy phases can be acquired via both translations 
and rotations of the local atomic orbitals through an inertial 
(laboratory) frame of reference: dynamics that may be ex- 
pected to be significant in free-standing graphene flakes. 

The Berry-phase mechanism involving orbital rotations, 
while less important in bulk semiconductors, due to the crys- 
talline rigidity, turns out to be of relevance in graphene. 
This two-dimensional (2D) material is embedded in a three- 
dimensional (3D) space, which allows for large out-of -plane 
displacements f9l (limited by surface tension and a weak 
bending rigidity) and thus strong dragging-like effects via 
rotations of the orbitals. Despite its single-layer structure, 
graphene is found to be one of the strongest existing materials 
IfTOl . with a high degree of elastic control due to its reduced di- 
mensionality. The electronic properties, in turn, are sensitive 
to the mechanical distortions fTTUMi, offering new pathways 
for controlling both the charge and valley degrees of freedom. 
In this Letter, we study coupling between the Dirac electrons 
and the time-dependent elastic distortions in graphene, for a 



general 3D excitation of the lattice. 

The crystallographic structure of graphene is that of a 2D 
honeycomb lattice of sp^ hybridized carbon atoms, with two 
inequivalent lattice sites per unit cell (triangular Bravais lattice 
with two basis points, A and E). The electronic properties of 
graphene are mainly due to the p, orbitals of the carbon atoms, 
which are perpendicular to the graphene sheet and give lise 
to the so-called n band. The other orbitals (sp^ hybridized) 
give rise to the cr bands, whose excitations are much higher 
in energy and the strong hybridization is responsible for the 
mechanical properties of graphene. The approximate tight- 
binding Hamiltonian describing the 2D graphene accounts for 
hopping between only the nearest-neighbor lattice sites and 
reads 10: 

Hg^YjYj [h^Rs)c^R,CR,^b, + H.c.j , (1) 

.! p=I,2,3 

where (cr) are the fermionic creation (annihilation) op- 
erators at site R, index s runs over all atom A sublattice 
sites, fp(-R) is the hopping amplitude to atom A at posi- 
tion R from its neighbor B atoms located at JS + bp, with 
bi = ( V3/2, l/2)fl, 62 = (- V3/2, l/2)fl, and ^3 = (0, -l)fl, a 
being the equilibrium interatomic distance and / running over 
all A sites. For the undistorted graphene, fp(-Rj) = fo, i.e., all 
tunneling matrix elements are equal, which allows one to find 
the spectrum easily by writing Hamiltonian ([T]) in the recip- 
rocal space: cjt a - Y^k ^k.o exp (ik ■ R)/ y/N, where a labels 
sublattice A or B, and is the total number of unit cells. The 
reciprocal space of graphene has hexagonal Brillouin zone, 
with two inequivalent points K and K' at its vertices. The 
resultant low-energy spectrum consist of two Dirac cones lo- 
cated at the K and K' points, around which the spectrum is lin- 
ear in momentum k, E(k) = ±Vf |A;|, with vf = 3toa^/2 being 
the Fermi velocity |8|. Transforming back to the real space 
description, the effective Hamiltonian around one of those 
points, say K point at (Att/ (3 V3a), O), is given by (putting 
^ = 1) //p = -ivfcr ■ V, where cr is a vector of Pauli ma- 
trices acting in the (A, B) sublattice basis. Since the K' point 
is related to K point by time reversal |5 1, their spectra, which 
define two valleys, are essentially identical. This Hamiltonian 
is responsible for many of the exotic electronic properties of 



2 



graphene, such as Klein tunneling fTSl, the half-integer quan- 
tum Hall effect 12j ,3J, zitterbewegung 1 16] etc. 

As a starting point, we generalize the above Dirac Hamil- 
tonian to nonuniform (but still real-valued) hoppings tp{Rj) + 
tQ. When for any position R,, the hopping tp(Rs) = tp, but 
tp + tpf for p p',we obtain the following Hamiltonian [9J: 

H'^ = vpcr ■ i-iV + A) , (2) 

withA^ = (l/3a)(fiH-f2-2f3)/foandAv = (1/ V3fl)(f2 - fi)/fo 
being the fictitious gauge potentials emerging as a conse- 
quence of the anisotropy in the hopping parameters tp. Ac- 
cording to the time-reversal invariance, which is unaffected 
by a static deformation, the gauge fields have opposite signs 
at the K' point. Making these gauge fields space-dependent, 
the Dirac electrons are subjected to a fictitious magnetic field 
Bk,k' - ± V X A, with opposite signs in the two valleys. 

There are two main physical mechanisms that lead to mod- 
ifications in the hopping tp(Rs) in static graphene: due to 
changes in distance between the carbon atoms and due to re- 
orientations of the atomic orbitals in a 3D-deformed graphene 
sheet. In the first case, the change in the hopping parameter 
tp(Rs), written in the WKB spirit, reads lH: 

6tf\R,,) - fo - 1} * -/?fo [ap(R,)/a - l] , 

(3) 

where ap(Rs) is the modified interatomic distance between 
the atoms at Rs and Rs + bp. To the leading order, ap{Rs)/a- 
1 ~ b'pbpe^ij/a^ (here and henceforth summing over repeated 
indices /, j, k,... that run over the x, y coordinates), where 
ejj = eijiRs) is the strain tensor at the position R^, defined 
by Cij - {djUj + djUi)/2 + (dih)(djh), with u being the in-plane 
and h the out-of-plane displacements. This expression is in- 
dependent of the type of orbitals we are dealing with, be it s, 
p, d, etc. The /? parameter, however, depends on the orbital, 
and for graphene fi ~ 2 for the p. orbitals |5|. The resul- 
tant gauge field components read A^l^^ - /SeijKjkieki/a, where 
Kxxx = 1, Kxyy - Ky^y - Kyy^ = "1 ^rc nonzcro elements of 
the rank-3 trigonal tensor [9 | and e,j is the rank-2 Levi-Civita 
tensor whose nonzero elements are fi, = -eyx - 1. Written 
explicitly, the vector potential becomes (9] : 

A^\r) = £ [e,,(r) - e„,(r), -2e,,(r)] , (4) 

which for a nonzero magnetic field thus requires an inho- 
mogeneous strain. Note that we are defining the strain ten- 
sor and, subsequently, curvature characterisics of a graphene 
sheet, along with all gauge fields, with respect to its intrinsic 
flat coordinate system (i.e., in the graphene's frame of refer- 
ence). 

The second mechanism affecting hopping strength involves 
only the out-of-plane displacement h(r). As already men- 
tioned, such distortions cause rotations of the atomic p^ and 
the hybridized sp^ orbitals, which become intermixed, thus 
modifying the tunneling Hamiltonian |5|. To quantify the 
amount of orbital mixture, we introduce the normal to the sur- 
face and the unit vector connecting neighboring atoms at the 
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FIG. 1. (color online). The sketch of a two-atom (A and B) unit-cell 
basis of a dynamical graphene lattice. The p- orbitals (blue dumb- 
ells) instantaneously point along the directions Na b(t) normal to the 
graphene substrate (gray), which rotate at angular velocities Na.b(t)- 
Here, b(T) is the instantaneous distance between the two atoms, and 
Ua(B)(t) are their linear velocities. 



position Rs, respectively IfTTI : 

z-Vh(Rs) , bp + z6hp(Rs) 
Ns = , , d,,p = , , , (5) 

Vl + [V/;(i?,)]2 V«' + [Shp(Rs)? 

where dhp(Rs) = h(Rs + bp) - h(Rs). In the flat graphene, 
we can write ifTSl to = {pzA^ ~ Eo\Pz,b) ~ -2.8 eV and fg = 
{py,A\H - Eolpy^B) ~ 8.4 eV, where Rb ^ Ra + ^3, H is the 
full microscopic single-particle Hamiltonian in the lattice, and 
E(, is the energy corresponding to the p orbitals in an isolated 
carbon atom. In a curved graphene sheet, this is generalized 
to tpiRs) = (AT, • Ps\H - Eo\Np ■ pp\ with p = (px,py,p,) 
being the p-orbital wave functions at a respective position in 
the lattice, i?, or Rp = R, + bp. We thus get for the change 
in tunneling due to curvature: 

6tf\Rs) = toiNs ■Np-l) + {t[, - to){Ns ■ d,,p)(Np ■ d,,p) 

^ - | [(bp ■ V)VhiR,,)f - [{bp ■ VfhiRs)f . (6) 

In this case too the changes in the tunneling tp lead to a fic- 
titious gauge potential A^'^' = lYjp^ijbpStf^ /3ti)a^. Written 
exphcitly, the gauge field components AfHr) read: 

Af^(r) = ^eij ^Kju6,„„ + ^Kjk„,i^ adkmh(r)di„h(r) , (7) 

where t = (f^ - to)/2to and the rank-5 tensor Kji^rsn has the 
following nonzero elements: Kyyyyy - -5, Ky^xxx - ■ ■ ■ = 
Kxxxxy = 3, snAKxxyyy = Kxyxyy = ■ ■ ■ = Kyyyxx = 1. Expanding 
this out, say, for the A^"^' component: 

- ^-a [{d]hf - {dih){d]h) - 2{dxyhf] . (8) 
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We identify the first term, oc (1 + t), with the effective gauge 
field derived in Ref. [VIj while the second term, oc t, is a new 
contribution that was previously overlooked in the literature 
(similarly, we obtain 2 terms forA^f'^) lfT9l . 

One could try to extend the above static theory to the dy- 
namical case by simply assuming the strain tensor e(r) and 
the height h(r) to be time-dependent in the above expressions. 
The dynamic distortions then give rise to a fictitious electric 
field E(t) = -djAir) (r stands for real time), in addition to 
the aforementioned magnetic fields B = 'VxA. The effect of 
such fictitious electric fields |20| have been studied recently 
in connection with both electron energy relaxation in carbon 
nanotubes [22] and electron charge pumping in graphene Il23l . 
This quasistatic picture alone would, however, be incomplete, 
as it neglects Berry-phase effects engendered by the motion 
of the atomic basis orbitals themselves in the dynamically dis- 
torted lattice. To capture such effects in graphene, we start by 
deriving the effective electron tunneling Hamiltonian i/eff be- 
tween two isolated A/B sites in a time-dependent framework. 

The full Hamiltonian of an electron in the potential 
y(r,T) created by the two atoms, A and B, is Hab{t) - 

12m + y(r, r), with m being the bare electron mass. 
Typically, one then solves the Schrodinger equation ap- 
proximatively by choosing a trial wave function |i/'(t)) - 
[ca(t)|(Aa) + cbWIiAb)] exp (-/^or), where \^a.b) = \Na,b ■ 
Pa.b) and ca,b respectively stand for the atomic wave func- 
tions and their amplitudes at time t, and derive an effective 
Hamiltonian H^s in the A, B subspace. However, these single- 
atom basis wave functions are not instantaneous eigenstates of 
the two-atom Hamiltonian Hab, which has nonzero matrix el- 
ements between the states \4'a,b) and the other orbitals outside 
this effective subspace. In the static situation, one can neglect 
such couphngs, as their effect on i/eff appears in second order 
in perturbation theory with respect to the inverse atomic en- 
ergy splittings, while, on the other hand, the direct coupling in 
the \il/A,B) subspace is first order and does not suffer from the 
energy mismatch. For the time-dependent case, in contrast, 
the intra-atomic electronic structure is itself perturbed by dy- 
namics, which requires one to systematically account for the 
associated excursions outside of the effective subspace \^a,b)- 

To capture such higher-order effects, we perform a more 
general time-dependent expansion: W{t)) - \\Ij{t)) + 
Tja=A,BMp„ c„,a{T)\na) cxp (-iE„T), with c„,„(t) , |n„>, and E„ 
being respectively the amplitude, wave function, and energy 
of electron in an excited orbital n. Substituting this expres- 
sion into the Schrodinger equation, /5t|(/''(t)) = HabW(j))^ 
leads to the following set of differential equations for c^ aij)'- 

ic„,a(j) = Yj c«'.-'(^)<««l'^'ifi - ^« - idr\K,)e'^^"-^"'^' . (9) 

n' ,ly' 

We will solve these equations only approximatively, in pertur- 
bation theory, assuming the time dependence of Hab is slow 
on the time scale ~ AE, with AE being the typical or- 
bital level splitting, so that the system stays mostly in a state 
spanned by the p orbitals. Then, taking c„,q.(t) » for excited 
orbitals, and |c^(t)P + |cb(t)|^ « 1, we obtain in leading order 
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{na\HAB - E„ - idAngi) ^nE„-Eo)T 

En - Eq 



F(t) , 
(10) 



where F{t) is a function that varies slowly on the time 
scale Tj. By substituting this expression in the equations 
for ca,b('!") and neglecting all the fast oscillating terms oc 
exp {i{EQ - E„)t\, we end up with the effective Schrodinger 
equation /c(t) - Hf-gcij), where i/gft 
6t = 5^^^ + 5f-"^ and 



f£v + 511.,., with 



6f ^\T) ^mt(T)U(T) ■ b(T) , (1 1) 

5f^"^(T) =^ {{Na ■ Pa\Nb ■ Pb) - {Na ■ Pa\Nb ■ Pb)) ■ 

(12) 

Here, U{t) = u(t) + zh{T) is the sum of in-plane and out-of- 
plane displacements of the center-of-mass of the two atoms, 
b(T) is the distance vector connecting the two atoms, f(r) = 
fo + Stir) K {iffA\p^/2m + V{r,T) - E()\\j/B) is the quasistatic 
tunneling amplitude discussed above [i.e., 6t - 6t^^ + St^''^], 
c{t) = (ca(t), cb(t)), and S are Pauli matrices that act in 
this A, B basis. Equation ( fTT) can be recognized as a local 
Galilean boost, for the derivation of which it is crucial to take 
into account corrections stemming from higher-orbital excur- 
sions, Eq. ( [Tol l |24|. Eq. ( fTS] ) is due to rotations of the lo- 
cal atomic p orbitals (see Fig. [T}: It emerges in the time- 
dependent picture as a Berry-phase correction to the tunneling 
term defined in Eq. (|6|. Note that the Galilean boost term can 
be re-exponentiated to give f — > f' = fexp(-/mf7 ■ b), as is 
expected according to the Peierls substitution in the moving 
frame of reference. 

We are now equipped to construct the full graphene Hamil- 
tonian in the dynamical case. The new (dynamic) corrections 
to the tunneling, are purely imaginary, while the qua- 

sistatic ones, 6t'^'''\ are purely real, which allows us to write 
ftot = tQ+St+i6t, with 5t = dt'^Usf"'^ and St = 6f^^+6^''\ The 
resultant effective Hamiltonian around the K point still has 
the form of Eq (j2]i, but with the total gauge field Aiot{r, r) = 
Air, T)+A{r, r) being the sum of the quasistatic and dynamic 
contributions, respectively, where A,- = 2Yjp £ijbp6tp/3tQa^ 
was evaluated already and A; = 2 2„ b' 6tp/3toa^ consists of: 



ill + oh dih + —Kijkh djkh 



(13) 



a 

2f^ 



S()5ij6ki + 



s,, - so 



-Kijki 



SoKijkdlm H Kijikn, 



dk,„hdjih, (14) 



where = 6ij6ki + dikSji + Sudjk, sq = {Pz.a\Pz,b} and 

'^0 - {PyAPy.B) (taking their values to be at 0.1 and -0.3, re- 
spectively) are orbital overlaps between the nearest-neighbor 
atoms for Rb - Ra + b^,. We find that in the presence of 
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FIG. 2. The radial projection of tlie fictitious gauge field, A,.(r) = 
A(r) ■ f, stemming from different mechanisms discussed in the main 
text, for a Gaussian bump of radius L = 100 nm, height hg = 10 
nm and u) = 10' s"'. From left to right: a) static strain-dependent 
component Af\ b) static curvature-dependent component A'f\ c) 
dynamic Galilean component Af'\ d) dynamic curvature-dependent 
component A*"'. All the gauge fields are expressed in units of cr', 
where a = 1 .42 A. 



both in-plane and out-of-plane time-dependent distortions, the 
effective Hamiltonian for graphene is complemented by the 
new, purely dynamic, gauge fields A^^^ and A^''\ which form 
the main result of our paper The leading-order terms oc h{r, r) 
in A*"^'"' appear formally as the dynamic analogue of the static 
teruis^ respectively, by substituting one partial spatial 
derivative djh — » cJj/i, with / = x,y, in the static fields. The 
next-order terms are proportional to the (time) derivative of 
curvature, oc dijh, which do not appear in our static expan- 
sion. Repeating the analysis for the K' point (either micro- 
scopically or by time reversal of the /T-point solution) shows 
that we need to substitute A — > -A, while A has the same 
sign and magnitude as at the K point. This has profound im- 
plications on the electron dynamics, as the gauge fields acting 
at the K and K' valleys are no longer equal in magnitude, and 
this asymmetry can lead to time-reversal-breaking observable 
electromagnetic properties, like charge currents and (quan- 
tum) Hall effect. 

As a particular example, we assume a time-dependent gaus- 
sian bump, h{r, r) = /iq exp {-r^ jL?) cos (w r), with /iq and L 
being the amplitude and the extension of the bump, respec- 
tively. In Fig. |2j we plot the resultant gauge field component 
Ax for the four different mechanisms: strain, Eq. (j4|, static 



curvature, Eq. (|7]i, Galilean boost, Eq. ( [T3| ), and dynamic cur- 
vature, Eq. ([14]). We see that for a given Dirac point {K point 
in Fig. [ij, the static gauge fields Af'"\r) are invariant under 
C3 rotations, and under C(y rotations complemented by overall 
sign change (which corresponds to time reversal), while the 
dynamic terms aJ.'^''^* are simply invariant under Ce rotations. 
Let us now consider the quantitative difference between the 
four different gauge fields (restoring K): 

|fo|L3 aML\L} ' 

where Ei - fp- jmL^ is the free-electron energy with a wave 
number corresponding to the distortion size L. While it may 
appear that the largest contribution is that arising from strain 
(oc yS), it is in fact overestimated here, since graphene is elas- 
tically extremely stiff for strain deformations ifTOl . Bending, 
on the other hand, costs much less energy ||5|, so that the low- 
energy height perturbations should favor profiles with mini- 
mal strain Il25l . While we treat the planar distortions [setting 
them to zero in Eq. ( [T5] l] and the height profile h{r, t) as in- 
dependent parameters, in practice, a self-consistent treatment 
of elastic properties would be necessary, by balancing bend- 
ing, shear, and compression energies 11711 . While addressing 
this problem in detail is important in view of quantifying the 
strength of the emergent gauge fields, such a calculation is 
beyond our immediate goals and will be presented elsewhere 
||26| . It suffices to say that, due to different spatiotemporal 
scaling (once the strain and bending are relaxed for a given 
height profile), we can easily envision geometric and dynamic 
limits where different terms in Eq. ([TSj would provide the 
dominant contribution. 
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